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Abstract 

We propose TBA integral equations for 1-particle states in the 0(n) non- 
linear (7-model for even n. The equations are conjectured on the basis of the 
analytic properties of the large volume asymptotics of the problem, which is 
explicitly constructed starting from Liischer's asymptotic formula. For small 
volumes the mass gap values computed numerically from the TBA equations 
agree very well with results of three-loop perturbation theory calculations, 
providing support for the validity of the proposed TBA system. 



1 Introduction and Summary 



The two-dimensional integrable 0{n) non-linear a-models have been the subject 
of a huge amount of study because they often arise in experimentally-realizable 
condensed-matter systems, because their integrability offers powerful theoretical 
methods applicable and because they share a lot of common properties with QCD, 
such as renormalizability, asymptotic freedom, dynamical mass generation, and in- 
stanton solutions. A particular interest of these models appears in lattice field 
theory, where because of their low dimensionality, the 0(n) non-linear cr-models are 
used to test new ideas. 

Many integrable quantum field theories (including the 0(n) non-linear cr-models) 
have been exactly solved "on-shell" i.e. their mass spectra and their factorizable 
S-matrices have been determined pp. Important progress was made towards the 
"off-shell" solution of the integrable quantum field theories by the so called Ther- 
modynamic Bethe Ansatz (TBA) technique [2; which proved to be very successful jS]. 
This approach allows one to calculate the ground state energy of an integrable quan- 
tum field theory enclosed in a finite box with periodic boundary conditions, provided 
that the "on-shell" solution is known. As the result of this method one obtains a set 
of coupled nonlinear integral equations (TBA equations), by the solutions of which 
the ground state energy can be determined. The knowledge of the finite size depen- 
dence of the ground state energy makes it possible to probe the scale dependence in 
integrable quantum field theories. 

The derivation of the TBA equations is cumbersome in the case of the 0{n) 
non-linear cx-models because of the difficulties appearing in the classification of the 
solutions of the Bethe Ansatz equations in the thermodynamical limit. This problem 
has been solved by formulating the 0(n) non-linear a- models as limits of certain 
perturbed conformal field theories [U IHl 13 IE] , where this classification can be done 
making possible the derivation of TBA equations for these models. The cx-model 
TBA equations are the limits of those of the corresponding perturbed conformal 
field theories. 

At the same time it seems to be interesting to calculate the full finite size spec- 
trum of the 0(n) non-linear cr-models, because these energies would provide inter- 
polation between the small volume (perturbative) spectrum and the large volume 
(massive) spectrum of the theory. 

To be able to achieve this, one needs to generalize the TBA technique for excited 
states. Although a lot of different methods have been worked out to obtain excited 
state TBA equations in different models CHI IHl 112 CSl El , a general method has 
not been discovered yet. The derivation of functional equations for the eigenvalues of 
certain mutually commuting transfer matrices is one of the most powerful method 
to attack the excited state problem in integrable quantum field theories [T3j . 
Knowing the analytic properties of these eigenvalues the functional equations can 
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be transformed into integral equations J3] which generahze the standard ground 
state TBA equations to excited states. This has been worked out for some simple 
perturbed conformal field theories directly in the continuum |11J or starting from 
the integrable lattice regularization of the model |14j . 

Unfortunately in the case of the 0(n) non-linear cx-models it is not obvious how 
to derive such functional equations (if they exist), and how to get information about 
the analytic properties of the functions in the case of excited states. 

As a partial solution of the problem, in ^Hj a simple method was worked out 
for proposing one-particle excited state TBA equations in the special cases of the 
0(3) and 0(4) non-linear cx-models. The purpose of our work here is to extend this 
method to the case of general 0{n). 

In this paper we propose an (infinite) set of TBA integral equations describing 
the 1-particle excitation energies in the 0(n) non-linear cr-model for even n, defined 
in a finite periodic 1-dimensional box. The equations are not derived from first 
principles and their validity is based on two major assumptions. First, we assume 
that the excited states are associated with the same Y-system equations as the ones 
proposed earlier for the ground state of the system |H|. Second, we assume that the 
analytic structure of the Y-system functions is independent (at least qualitatively) 
of the system size and can be read off from the infinite volume solution, which can 
be constructed using Liischer's asymptotic formula [1^. The first assumption can 
be proved for the sine-Gordon model [IH] and (at least for a subset of the physical 
states of the system) also for other simple models including the 0(4) cx-model. The 
second assumption is difficult to prove, even for the sine-Gordon model. However, in 
our previous studies [HI EE] we verified the validity of these assumptions by showing 
that the TBA equations obtained using them give finite volume energy values that 
agree very well with the results of Monte Carlo simulations and (for small volumes) 
with the results of perturbative calculations. In the same spirit, we study the finite 
volume mass gap values of the 0{n) models (for n = 6 and 8) numerically here and 
show that again, for small volumes, the agreement with the results of perturbation 
theory is excellent. This can be taken as an indirect proof of our assumptions. 
Indeed, since our starting point is the large volume asymptotic formula, the fact 
that the mass gap is correctly given by our equations also for small volumes is a 
very strong indication for the overall correctness of our approach. 

The paper is organized as follows. In Section 2 we recall the TBA integral 
equations describing the finite volume physics of the ground state of the model 
and write down the corresponding TBA Y-system equations. In Section 3, using 
Liischer's formula and assuming that the same Y-system describes the excited states, 
we formulate the 1-particle TBA problem. We find the solution in the infinite 
volume limit explicitly for 0(4) and 0(6) and conjecture the explicit form of the 
infinite volume solution for general rank. In Section 4 we prove this conjecture. 
In Section 5, assuming that the analytic structure of the solution is always the 



2 



same as in the infinite volume limit, we transform the 1-particle Y-system into TBA 
integral equations. The numerical solution of these equations is briefly discussed 
in Section 6 and the mass gap values are compared to perturbative results. The 
discussion of some technical details are relegated into the Appendices. Appendix A 
discusses a TBA lemma that enables one to translate TBA integral equations into 
TBA Y-systems and vice versa. Some useful formulae related to the 0(n) S-matrix 
are collected in Appendix B. In Appendix C the solution of the constant Y-system is 
given together with that of the related Q-system. Finally in Appendix D the details 
of the proof of Section 4 can be found. 



2 0(2r) TBA system 

In this section we summarize the TBA integral equations proposed by Fendley [Hj 
for the ground state energy of the 0(n) non-linear cx-models for even n, n = 2r. It is 
formulated in terms of a finite TBA system describing a certain perturbed conformal 
field theory model that becomes the 0(2r) a- model in the limit of infinitely many 
TBA components. 

The finite problem is characterized by the real pseudo-densities eo(6') and eni^{6). 
Here the upper index a = 1, . . . , r labels the nodes of the Dynkin-diagram (r > 2) 
and the lower index runs in the range m = 1, . . . , k — 1, where k > 2 is the Kac- 
Moody level characterizing the conformal limit of the finite system. The cr-model 
corresponds to the limit — >• oo j^. 

We now make the following definitions. 

Yo{9) = Lo = ln[l + Yo{9)] (2.1) 

and 

There is a symmetry relation 

et'\e) = e'^\e), m = l,...,fc-l, (2.3) 

which follows from the reflection symmetry of the Dr Dynkin-diagram. The ground 
state energy is given by 

Eo{R) = -^ I decosh^Lo(^), (2.4) 



L 



1 



1 + 



Y}^\e) 



a 
m 



k-l. 



(2.2) 
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where R is the size of the (one-dimensional) periodic box and M is the particle mass 
in infinite volume. We will also use the notation i = MR for the dimensionless size 
of the system. The pseudo-densities have to be determined from the TBA integral 
equations 



eo(^) =Mi?cosh^-^ r dun(^,^^^-^^^;e-u) l[^\u) 

a=l \^ 9 J 



(2.5) 



and for a = l,...,r, m = 1, . . . , k — 1 



+^S:Ii(^)+^S:Lm-E46aW( 



b=l 



(2.6) 



where the kernel function 7l{a,j3; 6) is given by ()A.8|) . g = 2(r — 1) is the Coxeter 
number of the Dr Lie-algebra and (by definition) 

Li^\9) = Li^\9) = 0, a = l,...,r. (2.7) 

Finally lab is the incidence matrix of the Dynkin-diagram, which is unity if the 
nodes a, b are connected and is zero otherwise. 

We can now use the TBA lemma discussed in Appendix A to transform the 
integral equations into TBA Y-systems. We first do this for the equations ()2.fij) . for 
which it is sufficient to use the special case of the lemma. We get for a = 1, . . . , r 
and m = 1, . . . , k — 1 

■{i + Yi:Ue)}{i + YiU0)}U 



with boundary condition 

Y(^\9) = Yt\9) = 0, a = l,...,r. (2.9) 
For ()2.5|) we have to use the lemma in its general form. We assume that 

l + r/"Ve):ANZ llmgl < ^ ~^)^ ^ a = l,...,r-2, (2.10) 

9 
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where ANZ indicates that the function is analytic and non-vanishing in the given 
strip. If this is satisfied then we can use the TBA lemma to transform ()2.5|) into 

a=l ^ ^ 

(2.11) 

The TBA Y-system fl2.8|) . ()2.1H) is illustrated in Figure 1. It belongs to the class of 
TBA Y-systems with one massive node {Yq in this case). 




Figure 1: Infinite Y-system corresponding to the 0(2r) non-linear cr-model 



3 1-particle Y-system 

In this section we conjecture the structure of the 1-particle Y-system, generalizing 
previous examples fHl in particular the case of the 0(4) non-linear cr-model 
(corresponding to r = 2). Our starting point is Liischer's asymptotic formula [T7] . 
For large i 

M 1"°° 

Eh{R)-HM / d^ cosh^e-^^°'^^^?/(e) (if = 0,1). (3.1) 

Here H = and H = 1 correspond to the ground state energy and first excited state 
energy of the 0(n) cr-model, respectively and for 

H = 0: y{9)=n, (3.2) 

n 

H = 1: y{9) = q{9 + '-^), where q{9) =Y, S^. (3.3) 



6=1 
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For H = ()3.H) is giving the first (universal) term of tlie virial expansion [T^ . 
wliicli is proportional to the number of particles (n), while for H = 1 the coefficient 
function y{6) can be obtained from the exact S-matrix of the model 

On the other hand for a TBA system with a single massive node Y, the analogous 
expressions are 

Eh{R)^HM / d^ cosh^ ln{l + i;(^)} (if = 0, 1). (3.4) 

271- J-oo 

Comparing this with Liischer's formula we conclude that asymptotically 

F.(^) ^e-^^°^^^l/(0). (3.5) 
In our case we conjecture that the 1-particle energy of the system is given by 

Ei{R) = M / d^ cosh^ ln{l + Fo(^)}, (3.6) 

where the 1-particle Yq satisfies the asymptotic condition 

ro(^) ~e-^™^'^^?/(^) (3.7) 

with the 1-particle y{6) defined in ()3.3|) and together with the other Y-functions Ym^ 
satisfy the (A; — > oo limit of the) Y-system equations ()2.8|) and ()2.11|) . Actually the 
structure of the equations becomes more transparent after introducing the variables 
Xo, Xm'* with rescaled arguments: 

lo(^)=Xo(f), ll'^H^) = Xi"^ (^) ■ (3.8) 

In terms of these fl2.8|) is rewritten as 

Xl^\x + l)Xl^\x -l) = {l+ 5„l(5al + 5a25r2)Xo{x)] 



{^ + xtUx)]{i + xtU^)]\[ 



X^ix) V (3.9) 



b=i 



1 + Xl\x 



for a = 1, . . . , r, m = 1, 2, . . . and ()2.1Hl becomes 

r-l 



Xo(x+z(r-l))Xo(a;-i(r-l)) = JJ |l+x{"^(x+i(r-l-a)) ||l+x("^(x-z(r-l-a)) |. 

a=l 

(3.10) 

In order to be able to transform ()3.9|1 and ()3.1()|1 into integral equations we need 
to know their analytic structure (position of poles, zeroes) in the strip relevant for 
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the application of the TBA lemma of Appendix A. Our strategy now is to study 
the large £ limit of the system first and try to establish the analytic structure of the 
Y-functions in this limit. We assume that in this limit 

Xo{x) ^ e — ), (3.11) 

9 

X^^^x) ^ Xt\x), (3.12) 
where Xm^ (x) is an ^-independent function. It follows that 

X,{x + ^ir-l))Xo{x-^{r-l)) ^(y + f - f ) (3-13) 

and using the results in Appendix B that 

r-l 

Jjj 1 + Xl""^ (a; + z(r - 1 - a)) } 1 1 + ^ (x - z(r - 1 - 

a=l 

4r^7(x + i{r — l))7(x — i{r — 1)) 
" {x^ + 4:){x^ + g^) ' 

where 

7(x)=x^+ (3.15) 
r 

The £ — > oo limit of ()3.9|1 simplifies to (a = 1, . . . , r m = 1, 2, . . . ) 



(3.14) 



X^\x+z)X^\x-z) = {l+xi:Ux)}{l+X^:U^)} n ( , ) , (3.16) 

1 + X^'(x) 



b=l 

which is the k oo limit of the standard o Y-system 

For the simplest case, r = 2, we have to determine the functions Xm^ for m 



0, 1, ... . In this case ()3.14j) determines 

Af (^) - ^ (3.1T) 
and using ()3.1fij) recursively leads to the complete solution 

x^ + {m + ir ^ + X^[x)- ^2 + (^ + 1)2 • (3.18) 

We note that this is of the form 

Xi^\^) = Vm , 1 + X^r^\x) = (1 + Vm) , ^ , (3.19) 

X^ + l+ym X^+l + Vm 



where Um are the constant solution of the r = 2 Y-system (see Appendix C). 

In the next example we are going to study, the r = 3 case, there are two sets of 
functions, and Xm ■, m = 0, 1, . . . . The general solution of ()3.14|1 is 

^ + ^""^ - 4 x2(x2 + 9) ^^-^^^ 

1 + Xf\x) = ^^±li;(x + ^)^Pix-^), (3.21) 
3 7(x) 

where iplx) is an arbitrary real analytic function. For the simplest choice, iplx) = 1, 
this can be extended to a solution for all m = 0, 1, . . . 

_ m(m + 4) {x^ + bl^,){x^ + bl_,) 
"^^"^^ " 4 x2(x2 + (m + 2)2) ' ^^-^^^ 

^i^Hx) - "^V'^ ^l. ^ (3-23) 
3 + 0^ 

_ (m + 2)2 (x2 + (&^ + l)2)(^2^(&^-l)2) 
^ "^^^ " 4 x2(x2 + (m + 2)2) ' ^^-^^^ 



where 



(m + l)(m + 3) 



(3.26) 



Again, we can recognize the structure 

2 2 1 

X!^\x) = ^ ^, l + Xi\x) = (l + l/i^)) / ^ (3.27) 

X2 + 1 + y)^^ X2 + 1 + y)^^ 

where ?/m are the constant solution of the Y-system. 

Based on these suggestive examples, we now formulate our conjecture for general 
rank. For the infinite volume solution for general r the functions associated to the 
"fork" part of the Dynkin-diagram are given by 

Xt'\x) = X«(x) = yt'^- (3.28) 

1 + Xt'Hx) = 1 + X«(x) = (1 + yt'^) "^'^Vd ' ^3.29) 

where ym~^^ are the constant solution of the Dr Y-system (see Appendix C). Starting 
from here, the Y-system equations ()3.16p determine the remaining functions Xm\x) 
for a = 1, . . . , r — 2 and m = 1, 2, . . . . Thus this solution, if exists, is unique. This 
will be discussed in the next section. 
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4 Infinite volume solution 



In this section we give the complete solution of the infinite volume Y-system explic- 
itly extending the conjectures ()H.28IIH.2??|l . We start by introducing 

yt^ : a=l,...,r m = 0, 1, . . . , (4.1) 

pt^ : a = 0,l,...,r-l m = 0, 1, . . . , (4.2) 

where y^m are the solution of the constant Y-system described in Appendix C and 
p'm are some additional (real) parameters. We also introduce the squares 

^t^ = {pt^f a = 0,l,...,r-l, m = 0,l,.... (4.3) 

Motivated by the r = 3 solution ()3.22|) and ()3.24|) we extend ()3. 28113. 2^ for 
a = 1, . . . , r — 2, m = 0, 1, . . . by 



and 



(4.5) 



We are now going to show that it is possible to chose the parameters pm so 
that ()3.28|) . ()4.4|) satisfy the Y-system ()3.16|) . This, together with the consistency 
between fl4.4|) and ()4.5p . requires 

e^'^ = l + Z/r'^ } m = 0,l,..., (4.6) 



m 


-1) 


= 


■{^- 
m 


-2) 


= 1 + v^'^-i) 


■(0) 
m 




(m + r — 1)^ 


.(a) 
.0 




(r - 1 - af 



0,l,...,r-l (4.7) 



and 

2(i+i/L"^)(d"^ + i) = z/L^ndti + d"-i) + d"^'^ + (4-8) 

fora = l,...,r — 2, m = l,2, Using the "boundary conditions" ()4.6H4.7|) . ()4.8p 

alone is already sufficient to calculate all other . We have an overdetermined 
problem and the question is if there is a consistent solution (which is unique, if 
exists). 
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In addition we have to require that Liischer's constraints ()3.14|) are also satisfied. 
This gives 



r-l 



n (1 + l/f"^) = 2r, (4.10) 



a=l 

(I) (r-2)(r^-l) 

er = ^ — j. — -■ (4.11) 

Using the formulae of Appendix C we can see that ()4.10|) is indeed satisfied since 

n (1 + Z/!"^) = gS'^ = 2r. (4.12) 

a=l 

With the help of and the Q-system elements Qm^ introduced in Appendix C we 
now define new variables 



„(a) _ Q(a) Ma) 


a = 0, 1, . . . , r — 1 


m = 0, 1, . . . , 


(4.13) 


f){a) _ Q{a) 


a = 0,l,...,r-2 


m = 0, 1, . . . , 


(4.14) 




m = 0, 1, . . . . 




(4.15) 



In terms of these new variables ()4.8|) and ()4.9|) together with the Q-system equations 
(jdljl can be written as homogeneous quadratic equations for a = l,...,r — 2, 
m = 1, 2, . . . as follows. 

+Qr'^Qr'\ (4.16) 



m+l Vm— vZi + QT'^ C"'^ + QZ-'^ vr'K (4.17) 

yim ) — 'lm-1 'lm+1 ^ 'Im 'Im ■ \^-^°J 

Further simplification is obtained if we write all equations in terms of the variables 
km and km , where 



Qrn ^rn ^m\ 1 ^ 0, 1, . . . , T 1 

(a) ^ Q{a) J^ia) ( m = 0, 1, . . . 

Im ^m ^m I ' ' 



(4.19) 



Here 

a 

G(^) = l, G'(:) = n^- a = l,...,r-l (4.20) 
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and finally 



j=b V / 



(4.21) 



From the solution of the constant Q-system we can read off that 



1,(0) _ 1 r,(l) 



m + r 



- 1 



(4.22) 



1 



m = 0, 1, . . . 



r 



and we can prove that if the variables km and km are identified with the ones dis- 
cussed in Appendix D then the quadratic equations ()4.1fjll4.lR|l are satisfied. ()4.1fij) . 
flTTTjl and (prTR|l follow from (ffT9|l . fflTTT]) and (ffTTn|l respectively. We can show 
by using the results in Appendix D that all the remaining requirements ()4.6|) . ()4.7|) 
and ()4.11|) are also satisfied. 

To summarize, we have shown that the solution of the infinite volume Y-system 
equations for the 1-particle case is indeed of the form ()H.28|1 . ()H.29j) . ()4.4|1 and ()4.5|1 
and calculated the constant parameters of the solution explicitly. 

5 1-particle TBA integral equations 

Having determined the infinite volume solution of the 1-particle Y-system equations 
we now turn to the problem of transforming these equations into TBA integral 
equations. 

We start by writing down the integral equations for Xm{x) associated to 
In this case we can use the special case of the TBA lemma and we only need to 
know the positions of the poles and zeroes (if any) of the functions Xm{x) (for 
a = 1, . . . , r — 1, m = 1, 2, . . . ) in the complex x strip with |Imx| < 1. From the 
explicit solution in the previous section and noting the inequality ()D.20|) we conclude 
that 



More precisely, we have established ()5.1ll5.3p for the infinite volume solution only. 
However, since in all other examples studied previously we found that the distribu- 
tion of poles/zeroes does not depend on the volume we shall assume in the following 
that ()5.1II5.H|1 hold for all values of the volume i and that no other poles/zeroes 
occur in the strip |Ima;| < 1. 

This structure can be taken into account by writing 



X^-^\x) has 
Xt'\x) has 
x['-^\x) has 



zeroes at x = ±i. 



double zero at x = m = 1, 2, . . . , 
double pole at x = m = 1, 2, . . . , 



(5.1) 
(5.2) 
(5.3) 



Xo{x) = e 



_£cosh^g/3o(x) 



(5.4) 
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and 



(forr>4) X';^\x) 



where 



e^™(^) a=l,. 



t{9) = tanh 



r — 3 



7l9 



m = 1,2, ... 



(5.5) 



(5.6) 



The right hand side of ()3.9|1 is regular for real x and has a double zero at x = only 
for a = r — 2, m = 1. Let us introduce 



£(;-i)(x) = In {1 + t\x) e^^'""(")} 
it-'\x) = \n{T\x) + e^-"'^^^} 
(forr>4) £(^)(x) = ln{l + e^™'(^')} a = l,...,r-3 



m = 1, 2, 



and 



where 



i^^^\x) = 61 Lq{x) - 25,"_2 In (1 + cosh — ) 



a = 1, . . . , r — 1, 



Lo(x)=ln{l + e-^™^'^^e*(^)} 



In the TBA equations we will use the combinations 



r-l 



^ini.^) — ^m+l(^) + ^m-l(^) + ^ah [/^m"* (^) ( 

6=1 



defined for a = 1, . . . , r — 1 and m = 1,2,..., where 



r-2r-l 



r-2 



1, 



(5.7) 
(5.8) 

(5.9) 
(5.10) 

(5.11) 



otherwise lab is the same as lah- 

Putting everything together, we can write down the 1-particle TBA equations 



Tlx 1 

13^^^ (yx) = Smi S^_2 In cosh J 



cosh 



Tr(x—u) 



(5.12) 



for a = 1, . . . , r — 1 and m = 1,2,.... Note that terms proportional to (5^_2 here 

and in ()5.8p take care of the zeroes at a; = ±z for '^^ {x) and simultaneously for 
the factor r^(x) on the right hand side of the a = r — 2m = l Y-system equation. 
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For the TBA integral equation associated to Xq{x) we need the TBA lemma in 
its general form. The following observations are necessary to be able to transform 
()3.10|1 into an integral equation. 

Xq{x) has zeroes at x = ±ip^p (in the strip |Imx| < t — 1), (5.13) 



1 + x["'\x) has zeroes at x = ±i(^p^^^ ~ ' ^ 



poles at X = ±ip\ 



(a+l) 



Imx| < r — 1 — a (5-14) 
(a = l,...,r-2), 



1 + : regular for real x. (5.15) 

(fH?T!?|l follows from (pTT|l . (lR6|l and (ICTl) . (fHm|l and (IKTTH|) follow from the explicit 
solution ()4.5|1 and ()3.29|1 . respectively, taking into account the inequalities ()D.20|1 
and (fiT22|l . As before, we can prove this analytic structure for the infinite volume 
case only. We assume the same analytic structure holds (at least qualitatively) for 
all values of the volume parameter i. 

In the notation of Appendix A we have 

r^ = a-l+p[^\ a = l,...,r-2 (5.16) 

for the complex zeroes and similarly 

Pa = a + p'i^^\ a = l,...,r-3 (5.17) 

for the complex poles^. The product of the relevant i?-factors is 

B{x) = n r f ^ + r ( ^-^i<^-_YP'^'^ ] (5.18) 



a=l 



and similarly for the i?-factors: 



a=l 



The prefactor in ()A.7|1 becomes 

B(x) 



Bix) ^(xWl'\^(--ip'l' 



^p^i = and consequently 1 + x'\^ (a;) has a double pole at x = 0. 



(5.20) 
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^ 


TBA 


PT 


Liischer 


0.001 


0.401699(7) 


0.40170 




0.003 


0.454323(7) 


0.45432 




0.01 


0.530922(3) 


0.53091 




0.03 


0.628234(3) 


0.62821 




0.1 


0.7881710(6) 


0.78810 


0.3159 


0.3 


1.0318308(5) 


1.03165 


0.6409 


1 


1.58529672(4) 




1.4128 


2 


2.3259656(5) 




2.2838 



Table 1: Results for the finite volume mass gap z in the 0(6) non-linear cr-model. 

Remember this prefactor is needed for a solution of the equation that is analytic 
and non- vanishing in the strip |Imx| < r — 1. The poles explicitly present in ()5.2()|1 
are canceled by the corresponding zeroes coming from the analytic continuation of 
the integral in the exponent. But according to ()5.13p we are looking for a solution 
that has zeroes at exactly the same positions! Thus in the solution for Xq{x) this 
prefactor is absent and the correct 1-particle TBA integral equation is simply 

/"OO 

/3o(x) = ^/ dun{r -l,r -l-a]x-u)6^\u). (5.21) 

a=l 

To summarize, the 1-particle TBA integral equations for the 0(2r) non-linear 
cr-model are ()5.12|1 and 1)5.2111 . The finite volume 1-particle energy in this model is 
given by 

Ei(R) = M / dx cosh — Lo(x). (5.22) 

2^ J-oo 9 

6 Numerical results 

In this paper we proposed TBA integral equations to calculate the finite volume 
1-particle energy in the 0(?t,) non-linear cr-models for even n. This, together with 
the corresponding equations for the ground state [SI, enables us to compute the mass 
gap of the model in finite volume exactly. Because we do not know how to derive 
these equations from first principles, the proposed equations remain conjectures. 
For this reason, it is extremely important to compare the resulting mass gap with 
the results of other approaches to the problem. 

We have computed the finite volume 1-particle energies together with the corre- 
sponding ground state values for a number of C. values numerically. The results for 
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the dimensionless mass gap 

z = R[E,{R) - Eo{R)] (6.1) 

as function of the dimensionless size of the system, i, are given in Table 1 and 2 for 
the 0(6) and 0(8) models respectively. The corresponding values calculated from 
Liischer's asymptotic formula ()3.1|) and three-loop perturbation theory j2l] are also 
given in the tables. 

The numerical computation proceeds as usual by iteration (after the discretiza- 
tion of the numerical integrals), starting from the £ — ^ oo solution. An extra com- 
plication for the numerical calculation of the mass gap from our TBA equations is 
caused by the presence of infinitely many unknowns. We have used here the same 
method as was applied to the case of the 0(4) model, which is described in ^| in 
some detail. This is not repeated here. Also the PT formulae necessary to compute 
the three-loop perturbative mass gap can be found in [16^. 

Here we just mention that we checked that to achieve the precision we quote 
it was sufficient to replace the original equations in the numerical computation by 
those of a cutoff problem, which is obtained from the TBA equations by assuming 
that all functions /5o(x), f3m\x) are equal to their infinite volume limit for x > 100 
or m > 140. For sufficient convergence, we need about 10-12 thousand iterations in 
these computations, which require about 7-8 hours of CPU time on a fast PC for 
the 0(6) problem, and almost twice as much for 0(8). 

From Table 1 and 2 we see that our numerical results agree with Liischer's asymp- 
totic formula for large volumes and agree with the results of PT for small volumes. 
Since our considerations started with Liischer's formula, the large i asymptotics is 
built in, but the very nice four-digit agreement between TBA and PT for small i is 
non-trivial. This agreement makes us confident that our TBA equations are indeed 
the exact answer to the question of the finite volume mass gap problem in non-linear 
(T-models. 
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TBA 


PT 


Liischer 


0.001 


0.38677(2) 


0.3868 




0.003 


0.43870(2) 


0.4387 




0.01 


0.51468(2) 


0.5147 




0.03 


0.61179(2) 


0.6118 




0.1 


0.772630(1) 


0.7727 


0.2486 


0.3 


1.019830(1) 


1.0200 


0.5536 


1 


1.584599(1) 




1.3465 


2 


2.333139(1) 




2.2656 



Table 2: Results for the finite volume mass gap z in the 0(8) non-linear (T-model. 

A TBA Lemma 

In this appendix we formulate a simple lemma |22j that we use to transform TBA 
Y-systems into integral equations and vice versa. 

Let the parameters a and (3 satisfy < /? < a and let f{9) be a (real) meromor- 
phic function in the strip |Im^| < (3 with no poles or zeroes at the boundaries. Let 
us denote by 

{rj} the set of real zeroes, 

{pj} the set of real poles, 

± ir'^^ the set of complex zeroes in the strip (0 < < /5), 

{pij ± ip'^} the set of complex poles in the strip (0 < p'^ < [3). 



Now we define 
and 

B{e) = n 

m = n 



r 



a J \ a 



a / \ a 



(A.l) 



Let us also introduce 

f^ = a-/3 + r'^ and p„ = a - /3 + p'^. (A.2) 



(A.3) 



(A.4) 
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where t{6) = tanh(7r0/4). 

(A.5) 

The function 

m = ^m, (A.6) 

where uj is the sign of f{oo), is positive along the real axis and we can define 

= ||y exp 1 1" du n{a, (3;^-u) In f{u) | , (A.7) 



with the kernel function 



1 cosh cos I? 



7^(a, f3;e) = - ^. (A.8) 

^ ^ a cosh ^ + cos ^ ^ ^ 

a a 

Now the TBA lemma, which can be proved by explicit integration, states that 
the function F{6) can be analytically continued into the strip |Im^^| < a, it is real 
analytic and non-vanishing (ANZ) there and satisfies the functional equation 

F{e + ia) F{e - ia) = f{e + f{e - (A.9) 

A special (limiting) case is /3 = 0. In this case the lemma states that if (f){6) is 
real analytic in the neighbourhood of the real axis then the function defined by 

is real analytic in the strip |Im^| < 1 and satisfies 

g{e + i) + g{e-i) = m- (A.ii) 
B 0(n) S-matrix 

In this appendix we study the properties of the function 

n 
6=1 

which appears in Liischer's asymptotic formula. The matrix elements of the 0(n) 
S-matrix are given by P 

Sf^{e) = ai{9) 5"^ 5"^ + a2{0) d'^'d^^ + a^{0) rt^^ (B.2) 
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where 



2ni9 S^^\9) 
' m - 9 {n - 2)9 - 2m 
S^^\9) 



{n - 2)9 



-2m 



{n - 2)9 - 2m' 
S^^\9) 



[n - 2)9 - 2m 



and 



-2-IT- 



5(2) (0) = - exp <! 2^ / — sin(^c^) '"^^ + ^" 



uj ' ' 1 + e 
From this we get for the combination (jB.lj) 

S^^\9)h{9-f) 



q{9) 



m) [{n - 2)9 - 2m\' 



where 



TT 



h{9) = n{n - 2)9' + — {n + 2){n- 4). 



Sphtting the Fourier kernel in ()B.4|) into two parts we can write 



S^'\9) = -K{9)M{9), 



where 



K{9) = exp <J 2z ^ ^ sin(^cu) 



and 



M{9) = exp{-i j — sin(6'cj) 



e 2 -I- e" 



7VUJ ^ 

27r — 



cosh^ 



B.8jl can be exphcitly evaluated: 



Ki9) 



Using this in ()B.5|) we get 



yi9) = q{9 + 



9 -in {n- 2)9 -2m 
9 + m {n- 2)9 + 2m 



M{9 + f) h{9) 
{9+'^) [{n - 2)9 + f{n + 2) 
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M{6) cannot be explicitly evaluated in terms of elementary functions but for our 
purposes it is sufficient to notice that it is an analytic and non-vanishing (ANZ) 
function in the strip < Im < tt and satisfies 

N 6 + 2711 (n — 2)6 + iirn 

M(6)M(6 + i7r) = —- ) , B.12 

^ ^ ^ ^ 6-iTT {n-2)6-27ii' ^ ' 



which can be shown using the Fourier representation ()B.9|) . It then follows that the 
combination relevant for the TBA Y-system is a rational function and is given by 

yy^^jyy-^)- {6-+7.^)[{n- 2f6- + M • (^-^^^ 

Finally we will compute the function (jB.llll along the real axis. We get 

y(0 ^ ;«>-^^°«' (B.14) 

Ve + + ^J{n - 2Ye + f (" + 2)' 

where 

Mo(0 =exp^J — cos(^u;) tanh — |^e"~ + e 2(n-2) J j . (B.15) 

C Constant Y-system and Q-system 

In this appendix we summarize the results for the standard Q-system [23] • We 
first study the generic case r > 3 and define for finite k > 2 the variables Qm^ 
where the upper index a = 1, . . . , r labels the nodes of the Dynkin-diagram and 
the lower index takes values m = —1, 0, . . . , k, k + 1. They satisfy the Q-system 
equations 

{Qm"=Qtw:i,+fi{Q^^^y'\ (c.i) 

b=l 

where lab is the incidence matrix of the Dj. Dynkin diagram^ for a = 1 , . . . , r and 
m = 1, . . . , k — 1. By definition 

Q(") = Qi^^ = 1, qLI = git =0, a = 1, . . . , r. (C.2) 

To simplify some formulae we will also use the definition = 1 (m = 0, 1, . . . , fc). 
Due to the refiection symmetry of the Dr Dynkin-diagram we have Qm~^^ = Qm- 



^i.e. lab = 1 if the nodes a, 6 are connected and vanishes otherwise 



19 



From the solution of the Q-system we can construct 



I N\ 2 



for a = 1, . . . , r and m = 0,1, ... ,k and it is easy to show that they satisfy the 
constant version of the standard Dr Y-system: 

r / (b) \ ^''b 

te^)' = (1 + + ylti) n 7^ (C-4) 

for a = 1, . . . , r and m = 1, . . . , k — 1. 

The unique positive solution of the Q-system ()U.1|) is characterized by [22] 

p(m + r- l) ^j^^p(n^_ m = -1. 0. . . . . 1. 1 + 1. (C.5) 

where 

= sin I — I (C.6) 
\k + gj 

and g = 2(r — 1) is the Coxeter number of the Dr Lie-algebra. The other components 
Qt^ for a = 2, . . . ,r are given by more complicated expressions, which can be 
obtained recursively using the Q-system equations and the result ()(15jl . 

For the 0(2r) TBA system we need the k —>■ oo limit of the Q-system and the 
constant Y-system. The solution ()C.5|1 in this limit simplifies to 



Oa. ^ !!i±I^ ^j^' !!i±l, „._i,o,.... (c.7) 

r — 1 7 

i=i 

In the first nontrivial case, r = 3, the solution of the constant Y-system is 

/IN m(m + 4) m(m + 4) ,^ , 

£^ = ^^^, £^ = ^^ m = 0,l,.... (C.8) 

The Q-system relevant for the r = 2 case is 

(Q^f = Q,n+lQn.-l + 1 (C.9) 

for m = 0, 1, . . . , A;, with boundary conditions Qo = Qk = 1, Q-i = Qk+i = 0. The 
Y-system is 

= + + m = l,...,k-l, (C.IO) 
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which is solved by 

ym = Qm+lQm-l, I + y.m = Q^n^ 171 = 0, 1, . . . , k. (C.ll) 

In this case the complete solution is 

= (C.12) 

P{1) 

which reduces to 

Qm = rn+1 and ym = rn{m + 2) (C.13) 

in the /c — > oo limit. 

D Linearization of the Q-system equations 

In this appendix we study a system of linear recursion relations that are useful in 
finding the solution of the infinite volume TBA Y-system. As a "byproduct" of 
this construction we find how to linearize the recursion relations for the infinite 
Q-system. 

The variables in these recursion relations are 

kl:\kl^\ a = 0,l,...,r-l, (D.l) 
where r > 2. The "initial conditions" are 

e = l, A;« = ^^, ¥^^ = {m + r-lf, (D.2) 
9 

where g = 2{r — 1) , and the rest of the variables are determined recursively by the 
relations 

(2m + ^7)4") = {g- a)kt'''^ + ak^'^ (D.3) 

and 

r(r - 1)A;(^) = (r - a - l)(m + r - l){(r - a)(m + r - 1)A;^) - ak^^~^^}. (D.4) 



()D.3|) is valid for a = 1, . . . , r — 2, while ()D.4|) for a = 1, . . . , r — 1. It is clear 
from these relations that km^ and km^ are polynomials in m (of order a and a + 2 
respectively). 

Using the definition ()D.3|) we can prove by induction that the variables km 
satisfy further linear relations. These are 

{g - 2a)kt^ = {g + m)kt\ - mk^U (D.5) 
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for a = 0, 1, . . . , r — 1, 

^m + g- a)k[^} = {g + m)k^:l, + ak^-'^ (D.6) 
for a = 1, . . . , r — 1 and 

{a-g- m)kt) = mk^l, + {a- g)kt'-') (D.7) 

for a = 0,l,...,r — 2. Further these relations can be used to show that k^ can also 
be calculated from 

r{r — \)k^ = (r — a — l)(m + r-l){(r-a-2)(m + r-l)A;^) + (^-a)A;^+i)}. (D.8) 

Our variables also satisfy three sets of quadratic relations (for a = l,...,r — 2). 
These are 

^m + a){m + g- a){kt^f = m{m + + «(^? - a)k^^+'ht'\ (D.9) 



^rn + a){m + g-a){k^:,^ = m{m + g)kl:l,k^^^^^ (D.IO) 
and 

2(m + a){m + g- a)kt\kt^ + k^^) = m{m + g)[kt\,ktl, + ktl^k^U 
+ a{g — a)[k^'^^'^ k^^ + ^m'''^''^m ^'']- 

We can prove these relations by expressing all variables occurring in them in terms 
of fcj^^i and fc^^"^^ using the linear relations (ffT8|l . (frT4|l . (frT5|) . (flTfi^ and (fiT8|) . 

As functions of these two variables the quadratic relations ()D.9j) , (jD.lOjl and (jD.lljl 
are identities. 

In the main text we need the following special cases. 

= 1, fcj"^ = (r - a - 1)2 a = 0, 1, . . . , r - 1, (D.12) 

^ ^^^)=o, (D.13) 

JJ^) = {r-2){m + r){m + r-2) ^ ^^^^^^ 

kt'^ = -^{im + r-l)kt'^-im + g)kt_l^}, (D.15) 
r — 1 

e^) = (D.16) 

r — 1 



22 



For m > integer we define the ratio 



From ()D.12IID.'l6|l we now have 

r 

ej"^ = (r-a-l)^ a = 0,l,...,r-l. (D.19) 

The variables also satisfy the following inequalities. 

d")>l, a = 0,l,...,r-2, m>0, (D.20) 

where equality is for the case a = r — 2,m = only. 
Also 



(r — a — l)^(m + r — l)(m + r — 2) 



^ Sm ■ 



(r — a)(m + r — 1) 



2 



(D.21) 



r(r — 2) 

for a = l,...,r — 2,m>l. An important consequence of the above is 

r - a - 1 < < r - a, a = l,...,r-2. (D.22) 
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